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ABSTRACT 

Context. The Orion Nebula and its associated young stellar cluster are located at the front-side of the optically thick OMC-1 molecular 
cloud. In order to disentangle the cluster members from background contamination, it is important to know the extinction provided 
by the OMC-1, which is poorly known, the available measurements yielding contradictory results. 

Aims. Our main goal is to derive a new extinction map of the OMC-1, obtaining information about the structure of the OMC-1 and 
the Orion Nebula Cluster. 

Methods. The most recent near-infrared catalog of stars is used to study the distribution of reddening across a 0.3 deg 2 area covering 
the Orion Nebula Cluster. On the basis of the observed (H, H - K s ) diagram, we establish a criterion for disentangling contaminants 
from bona-fide cluster members. For contaminant stars, interstellar reddenings are estimated by comparison with a synthetic galactic 
model. A statistical analysis is then performed to consistently account for local extinction, reddening and star-counts analysis. 
Results. We derive the extinction map of the OMC-1 with angular resolution <5'. We also assemble a sample of candidate cluster 
members, for which we measure the extinction provided by the nebular environment. These extinction measurements are analyzed 
similarly to the contaminant sample, and an extinction map of the Orion Nebula is derived. 

Conclusions. The extinction provided by the OMC-1 is variable on spatial scales of a few arcminutes, while showing a general 
increase from the outskirts (A v ~ 6) to the direction of the Trapezium asterism (Ay > 30). The Orion Nebula extinction map is more 
irregular and optically thinner, with A v of the order of a few magnitudes. Both maps are consistent with the optical morphology, in 
particular the Dark Bay to the north-east of the Trapezium. Both maps also show the presence of a north-south high-density ridge, 
which confirms the filamentary structure of the Orion molecular complex inside which star formation is still taking place. 

Key words. ISM: clouds - Interstellar medium (ISM), nebulae - dust, extinction - open clusters and associations: individual: Orion 
Nebula Cluster - Methods: statistical - Techniques: photometric 
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OO 1. Introduction and to explore its spatial variations with the distance from the 

*n | cluster center, it is necessary to remove the contribution of non- 

OO - The closest event of massive star formation is occurring at cluster sources. In principle, this requires the acquisition of thou- 

O ' P resent in the direction of the Galactic anticenter (I = 209°, b = sanc ? of spectra of faint sources distributed over the bright neb- 

t-H i - 19°), in the Orion constellation. The Orion Nebula (ON, M42) ular background. On the other hand, the number of contaminant 

. represents the most spectacular signature of star formation ac- sources , both galactic and extragalactic, can be estimated using 

L; ■ tivity in this region. The ON is a blister HII region carved me most recent models for stellar and galaxy counts at various 

. ^ ! into the OMC-1 giant molecular cloud by the UV flux emit- wavelengths. The main complication in this case arises from the 

^ . ted by a handful of OB stars, the so-call ed Orion "Trapezium" pre sence of the OMC- 1 , which provides a backdrop to the ONC 

?H ; dMuench et alj |2008|; jO'Dell et alj|2008|). The Trapezium stars f high and non-uniform extinction. Deriving an accurate extinc- 

P_ ; are t he most massi v e mem bers of a rich (« » 2000 mem- tion map wouk j be beneficial not only to better discriminate the 

bersjMuench et afj d2002|)) cluster of young (1-3 Myr old, Q NC membership, but also to understand the 3-D distribution 
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|Da Rio et alj 12010)) Pre-Main-Sequence objects (Orion Nebula f t h e cluster, still partially embedded within the ONC, and the 

Cluster, ONC). Given its youth, vicinity, and low foreground ex- evolutionary history of the region, 
tinction, the observed luminosity function of the ONC can be 

converted int o a true initial mass fu nction with relatively modest i n this paper we present a reconstruction of the OMC- 

assumptions dMuench et alj|2002j and references therein). It is 1 extinction map based on the analysi s of th e recent near-IR 

largely for this reason that the Orion Nebula and its associated (NIR) source catalog of Ro bberto et alj d20Tob (R10 hereafter), 

cluster are regarded as a critical benchmark for our understand- r n Section |2] we briefly review the previous studies relevant to 

ing of the star formation process. me Orion Nebula region. In SectionSwe illustrate our statistical 

In order to build a reliable luminosity function, especially in method to disentangle background stars from the cluster pop- 

the substellar regime (brown dwarfs and planetary mass objects) ulation. By combining an estimate of the interstellar extinction 
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affecting each contaminant star with the source count density, we 
derive the OMC-1 extinction map. In Sect. |4] we apply a similar 
statistical procedure to the candidate cluster members, deriving 
an extinction map for the dust in the foreground Orion Nebula. 
Finally, in Sect. [5]we compare our maps to previous studies, and 
we briefly discuss their main features and similarities to stellar 
distributions. 

2. Overview of previous studies 

A number of previous studies provide results relevant to the issue 
of t he galactic red dening in the direction of OMC-1. 

ISchlegel et all (fl998l) (SFD98 hereafter) combined the 
COBE/DIKBE observations (100 //m and 240 yum) and the 
IRAS/ISSA observations (100 jum) to obtain a full-sky 100 /mi 
map with ~ 6' resolution. On the basis of the correlation be- 
tween the Mg line strength and the (B - V) color of elliptical 
galaxies, they were able to calibrate their column-density map to 
a E(B - V) color excess map. Their all-sky reddening map, avail- 
able through the NASA/IPAC Infrared Science Archival Dust 
Extinction Service web page, represents the benchmark for the 
following studies of our region. 

The accuracy of the SFD98 maps has been analyzed by 
lArce & Goodman! (l999), who compared the extinction map of 
the Taurus dark cloud with the extinction maps obtained using 
four other methods, i.e. 1) the color excess of background stars 
with known spectral types; 2) the ISSA 60 and 100 //m images; 
3) star counts; and 4) optical color excess analysis. These four 
methods give similar results in regions with Ay < 4. Their com- 
parison shows that the extinction map derived by SFD98 tends 
to overestimate the extinction by a factor of 1.3 -r 1.5. They as- 
cribe this discrepancy to the the calibration sample of ellipticals 
used by SFD98, which containing few sour ces with Ay > 0.5 
may l ead to a lower estimate of the opacity. lArce & Good man 
d 1999 !) also find that when the extinctions shows high gradient 
(> 10"' /deg, see their Fig. 1), its value is generally smaller than 
that given from SFD98, and argue that the effective angular res- 
olu tion of the SFD98 ma p could be somewhat larger than 6'. 

iDobashi et alj (2005) obtained another all-sky absolute ab- 
sorption map applying the traditional star-count technique to the 
"Digitized Sky Survey I" (DSSI) optical database, which pro- 
vides star densities in the range ~1— 30 arcmin~ 2 . Comparing the 
SFD98 maps to their results in the Taurus, Chameleon, Orion 
and Ophiucus complexes, they found that SFD98 extinction val- 
ues are generally up to 2-3 times larger (see their Figs. 46, 47 
and 48). A least-squares fit of the two extinction estimates in the 
A y < 4 gives a proportionality coefficient of 2.17, while the same 
analysis in the Ay > 4 suggests a coefficient of ~ 3. The authors 
discuss the possible causes of discrepancy. First, they stress that 
the SFD98 map are mostly sensitive to the total dust along the 
line of sight, while their maps measure the extinction due to the 
nearby dust, as the optical thickness is much larger in the visi- 
ble than in the far infrared. Second, they point out that the low 
resolution (~ 1 deg) of the temperature map adopted by SFD98 
cannot reproduce the typical high spatial frequency temperature 
variations across dark dense clouds. Finally, they suggest that 
SFD98 do not account for enhancement in the far infrared emis- 
sion by fluffy aggregates: this excess emission could explain the 
inconsistency between SFD98 map and maps derived using other 
methods. Unfortunately, DSSI optical plates saturate in proxim- 
ity of the Trapezium cluster, so a large fraction of the ONC field 
is excluded from their analysis. 

1 http://irsa.ipac.caltech.edu/applications/DUST/ 



An extinction map limite d to the inner ~5'x5' regio n of 
the ONC was provided by Hillenbrand & Carpenter (2000) on 
the ba sis of the C ls O column density data of iGoldsmit h et al.l 
(119971) . with a 50"spatial resolution. A comparison of the ex- 
tinction map they derive with the extinction obtained by SFD98 
shows that that the SFD98 extinction is generally ~ 3 times 
larger. 

In summary, the SFD98 extinction map for the OMC-1 is 
still the only one covering the entire ONC field. However, its 
spatial resolution (6') is limited and the accuracy, in a region as 
complex as the OMC-1, remains questionable. 

In Sect. [3] we overcome most of the issues listed above us- 
ing photometric data in the NIR bands, as they offer several well 
known advantages. First, 90% of the galactic population is made 
up of M dwarfs which, by virtue of their SED, radiate mostly in 
the NIR. Second, the extinction is lower at IR wavelengths and 
therefore stars here can be more easily detected through optically 
thick clouds. Third, the high surface density of field stars allows 
to overcome the limited spatial resolution typical of far-infrared 
wide-field surveys used to measure diffuse dust emission. As a 
caveat, however, we anticipate that in our the number of detected 
background stars will be limited by the brightness of the back- 
ground, wich affects the completeness limit of the NIR survey 
(Sect. [3). This effect will be taken into account in our analysis. 

3. The OMC-1 extinction map 

To compute a new extinction map of the OMC-1 region, 
we use the near-infrared photometry obtained by R10 at the 
CTIO/Blanco 4 m telescope in Cerro Tololo with the ISPI im- 
ager. The assembled photometric catalog contains JHKs pho- 
tometry in the 2MASS system for 7759 sources, 6630 labeled 
as point-like sources, spread over an area of ~30'x40' roughly 
centered on ^Ori-C (RA=05 /! 35 m 16.46 s , DEC=-05°23'23.2"). 
While this catalog is slight ly shallower than previous st ud- 
ies of the clus t er cor e (e.g. [Hillenbran d & Carpenter! d2000l) or 
iMuench et al.l (l2002l) \ it covers a sky area ~50 times larger, 
roughly centered on the Trapezium stars. We concentrate our 
analysis on the (H,H-Ks) color-magnitude diagram (CMD) of 
the ~6000 point-like sourc es with both H and Ks mea sured pho- 
tometry in the R10 catalog. Lomb ardi & Alvesl (12001) show that 
the inclusion of the / magnitudes can reduce the noise of the ex- 
tinction measurements by a factor of two, especially in regions 
characterized by low extinctions. Unfortunately, this is not our 
case. By using the full sample of stars detected in the H and Ks 
bands, we greatly increase the stellar density and therefore maxi- 
mize the angular resolution of our map, since the R10 photomet- 
ric catalog is shallower in the /-band. The R10 survey also lacks 
deep /-band observations of the North-East corner of their field, 
which makes the deep field covered by the the /-band smaller by 
about 9% than the field covered in H and Ks -bands. Later in this 
paper we will include the /-band data, when we will consider the 
extinction map toward the ONC, which lays in the foreground of 
the OMC-1. 

The {H,H-Ks ) CMD in Fig. Q]shows a characteristic bimodal 
distribution: a first group of stars (the ONC) is clustered at H- 
Ks-0.5 and H=12, whereas a second group of fainter objects 
appears clustered around H~17. This second group also has a 
peak, but this is a selection effect, as the number of faint sources 
drops at H ~Z 17.5 in correspondence of the sensitivity limit of 
the R10 survey. 

Figure Q] also shows a 2 Myr isochrone appropriate for the 
ONC (solid green line) and the density contour of the galac- 
tic stellar population along the line of sight of the ONC (solid 
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contours). The galactic pop ulation has bee n derived from the 
Besangon galactic model (Robin et al. 200310, computed at the 
galactic coordinates of the ONC over the RIO survey area of 
0.329 deg 2 . For each star we computed synthetic photometry 
in the 2MASS sys tem using the grid of atmosphere models of 
lAllard et"ail (12010). Since the galactic model also provides for 
each star the distance and interstellar extinction, both parame- 
ters were accounted for in the synthetic photometry. 

Of the 4868 sources with H<18 returned by the model, 4616 
(»95%) lie at distances larger than 420 pc. Their colors will 
generally be further reddened by the OMC-1. In Fig. Q] we plot 
(dashed lines) the galactic model reddened by Ay = 6, an ad-hoc 
value chosen to shift the contours in the same region occupied 
by the second peak at H~\l in the {H,H-Ks) diagram. As an- 
ticipated by RIO, this confirms that the second peak in Fig.[T]is 
fully compatible with the reddened galactic population. 

The peculiar morphology of the CMD, characterized by two 
well separated peaks associated to the ONC and the reddened 
galactic populations, allows us to assume that sources fainter 
than H-15 or redder than H-Ks = l3 (respectively below and to 
the right of the dotted lines in Fig. [TJ are largely background 
field star, being either too faint or too red to belong to the ONC, 
and vice versa for the other sources within the area. This is a 
first order criterion that can be further refined. Firstly, we can 
take advantage of the finer angular resolution of the ACS images 
obtained for the HST Treasury Program to discard 73 double 
stars and 166 extended obj ects, either background galaxies or 
proplyds ( Ricci et al.ll2008l) . Moreover, by looking at the JHK S 
color-color diagram for the subsample of stars with available / 
magnitude, we find 61 "field" stars well compatible with the lo- 
cus of reddened circumstellar disks (Section 3J. These sources 
have been therefore moved to the cluster sample. Our final sam- 
ples are thus made up of 1716 cluster sources and 1913 back- 
ground galactic stars. 

Figure [2] shows the projected density contours for the two 
samples of background and cluster stars. The relatively low 
surface density of background stars in the direction of the 
Trapezium cluster indicates that the optical depth of the OMC- 
1 reaches its maximum behind the Trapezium. Vice versa, the 
density of bona-fide cluster members reaches its maximum in 
the direction of the Trapezium cluster and decreases towards 
the outskirts of the Orion Nebul a, as already suggested by e.g. 
iHillenbrand & Hart mann (1998). The density of stars, either 
background stars or cluster members, is critical in our statistical 
derivation of the extinction map of the OMC-1 (Sect. l3~2l and of 
the Orion Nebula (Sect. |4| and drives the angular resolution of 
our derived maps. 

It is clear that the two samples are not pure, e.g. there are 
brown dwarfs and/or planetary objects belonging to the cluster 
that may have been improperly included in the background sam- 
ple. According to the previous stu dies of the ONC initial mass 
function (see iMuench et al.l 12008. and references therein), the 
fraction of sources having substellar masses decreases with de- 
creasing mass and therefore one may argue that the background 
sample is poorly affected by the presence of ONC member with 
very low masses. By constraining the stellar density to match the 
average extinction, we will further reduce the cases of spurious 
identifications. In any case, contamination will remain a source 
of uncertainty. 

While it is possible to assume that the intrinsic color of 
a cluster star (neglecting disk excess emission) roughly corre- 
sponds to that of a source on the 2 Myr isochrone, for the back- 

2 http://bison.obs-besancon.fr/modele/ 
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Fig. 1. Density map in the (H,H-K S ) CMD of all point-like 
sources in the RIO catalog with measured HK$ photometry 
(c olor-filled areas), together with the 2 Myr isochrone computed 
by IScandariato "et al. I din preparation!) (solid green line, the dia- 
mond symbols correspond to 1M and 0. 1M from top to bottom 
respectively). Two main populations are apparent: the ONC clus- 
tered at H=12 and the bona-fide background stars at H~17. The 
latter population is to b e compared with the synthetic galactic 
population provided by Robi n et ail (120031) (solid contours) after 
accounting for the extinction provided by the OMC-1. The case 
of a constant Ay=6 across the cloud is shown in dashed con- 
tours, the highest ones roughly overlapping with the observed 
background population. The dotted line represents our safe se- 
lection criterion: any source fainter than H-15 or redder than 
H-Ks = 13 is considered a field star, provided that its colors are 
not compatible with the reddened Meyer's locus (Sect. |4]i. 



ground population there is no such isochrone. In the following 
sections we describe the statistical approach (Sect. [3TTb we have 
used to derive the extinction affecting each background star and 
consistently derive the extinction map of the OMC-1. 

3.1. The extinction affecting each background star 

The basic idea of our method, similar to the one outlined by 
IHillenbrand & Carpenter! d2000l) . is illustrated in the left panel of 
Fig. [3] Each star is represented in the CMD by a 2D probability 
distribution (the ellipse at the bottom right of the figure), where 
the shape of the distribution is due to the correlation between 
the H and H-K$ measures, with their uncertainties. This density 
distribution has been uniformly binned in a 20x20 grid, resulting 
in 400 grid points. 

We project each grid point backward in the reddening di- 
rection, tracing a stripe which intersects the galactic population 
model. Each grid point is thus associated to a certain number 
of galactic stars with different extinction values, without any re- 
striction on the maximum extinction allowed. Each of these stars 
is a candidate to represent the dereddened photometry of the grid 
point. The central panel of Fig.[3]shows the cumulative distribu- 
tion F(Ay) of the reddening distribution of the candidates asso- 
ciated to the grid point shown in the left panel. 
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Fig. 2. Projected surface density of bona-fide background stars (left panel) and of bona-fide cluster members (right panel), as 
selected in Sect. [3] The contour labels give the density of stars in arcmirT 2 . In the direction of the Trapezium cluster (the star 
symbol indicates the position of ® l OriC, assumed to be the cluster center) we find the minimum density of background stars, while 
it increases towards the outskirts of the cloud. This is indicative of the decreasing optical thickness of the OMC-1 with increasing 
projected distance from the cluster center. Conversely, the surface density of cluster members steeply increases with decreasing 
distance from & l OriC. 



i 1 6.5 




15 



H-K s A v 

Fig. 3. Left panel - Example of the dereddening algorithm for the star with ID=578 in the RIO catalog. The star is represented in the 
CMD by a regularly gridded 2D gaussian density distribution, and to each grid element it corresponds a stripe of the galactic model 
(represented in dots) parallel to the reddening direction. Dereddened candidates for the outlined grid point are coded with colors 
from yellow to red according to increasing Ay. Central panel - Cumulative distribution F(Ay) for the colored synthetic galactic 
stars in the left panel. For each A y, this function represents the fraction of candidate dereddened photometries with extinction < Ay. 
Right panel - Weighted probability density function f(Ay) of the 40000 probabilistic Ay estimates for the representative case shown 
in the left panel, computed using a gaussian kernel. The weighted mean of the sample (dashed vertical line) is our best estimate of 
the true extinction affecting the star. 



Using the cumulative distribution of Ay, for each grid point 
we randomly draw 100 Ay-values and assign to the full set of 
100 estimates of Ay a weight given by the 2D probability den- 
sity distribution of the given grid point. Repeating this process 
for all grid points, we obtain a total of 400x100=40000 prob- 
ability estimates, each one with its own weight: the weighted 
mean of this sample gives the most likely estimate of the ex- 
tinction affecting a given observed star. The right panel in Fig. [3] 



shows the weighted probability function f(Ay) corresponding to 
the 40000 random estimates of Ay; to avoid data binning and 
for the sake of smoothness, the f unction is c omputed using a 
gaussian kernel density estimator ([Silverman! [19 86). The com- 
putation is repeated for all background sources, providing an ex- 
tinction estimate, with relative uncertainty, for each star of the 
sample. 
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3.2. Computation of the extinction map 

To compute the OMC- 1 extinction map, we project over the orig- 
inal survey area of RIO a grid of points every 100 x 100 pixels 
(corresponding to 30" x 30"). For each grid point we select the 
20 closest background stars within 1,000 pixels (5')- This limit 
is meant to prevent that, especially in regions with high extinc- 
tion (Ay >15) where the density of background stars can be 
extremely low, the algorithm collects sources over an area too 
widespread. Thus, by setting an upper limit to the distance of 
background stars, we tend to preserve the angular resolution of 
the extinction map at the price of working with smaller statis- 
tical samples and therefore with locally more uncertain results. 
The value of 5' has been tuned to make sure than even in the 
regions with the highest extinction and lowest background den- 
sity the algorithm will find a few background stars within that 
distance. This limit is also comparable to the scale size of major 
nebular structures seen in the NIR images. 

Each sample of 20 (or less) background stars is averaged us- 
ing an iterative sigma-clipping routine, deriving a robust esti- 
mate of both the mean local extinction Ay and standard deviation 
cr A - v . This would represent the solution to our problem, except 
that we want to make sure that the observed surface density and 
brightness distribution of the background sample is still compat- 
ible with the model. To this purpose, we apply the Ay extinction 
to the synthetic galactic model and compare the local density of 
the reddened model population to the observed one. 

The predicted density of the reddened galactic population 
has to be adjusted for the completeness of the survey, which 
statistically accounts for the success/failure of detection based 
on the magnitude of both stars and local background. In R10 
we computed the completeness levels in three concentric regions 
centered on 9 1 Ori-C, neglecting any variation inside each region. 
For this study we improve the completeness estimates using a 
locally-computed completeness level, as described in Appendix 
lAl for each grid point in the extinction map, our new approach 
allows us to compute the appropriate completeness adjustment. 

It turns out that, particularly within a few arcminutes from 
the Trapezium cluster, the number of predicted background stars 
is lower than the number of observed ones. The excess sources 
are, most probably, ONC members wrongly identified as back- 
ground stars, having photometry matching our background se- 
lection criteria. To statistically reject these extra sources, for 
each group of 20 stars we discard the background star with the 
lowest extinction Ay m ,„, assuming it is a cluster member, and re- 
place it with a new, further background star with Ay > Ay ml „, if 
existing (we are still constrained by the 1,000 pixels maximum 
distance). By decreasing the density of background sources, we 
obtain an extinction map which matches both the typical red- 
dening and source density of the galactic component, averaged 
over 5', or less. The result is shown in Fig. |H where we have 
interpolated the discrete gridding of the map using the IDL 
GRIDDATA.pro routine. The angular resolution ranges between 
1.5' and 5', depending on the density of background sources 
(Fig. [2). In correspondence of the highest extinction values the 
map reaches its worst resolution, as expected. 

The error map shown in Fig. [4] shows a typical uncertainty 
Ay <1. The regularity in this map is due to the fact that for each 
grid point, apart the Trapezium region, the algorithm collects the 
requested 20 stars within a distance shorter than the limit of 5'. 
It follows that the number of stars in the analyzed subsample 
(corresponding to the given grid point) does not depend on the 
stellar density, but on the combination of the density and the 
angular resolution of the extinction map. 



3.3. Biases in the selection of background stars 

3.3.1. Foreground contamination 

In deriving our extinction map we have assumed that the fore- 
ground contamination is closely modeled by the Besangon syn- 
thetic galactic population. It is otherwise well established that 
the ONC (also known as Orion subgroup Id) is neighbored 
by three somewh at older subgroups of stars (la, lb, and Ic, 
iBrown et al.ll 19941) . which are located at distances ranging from 
~360 to 400 pc and have ages from ~2 to 11 .5 Myr. The most 
likely subgroup from which we would see contamination in our 
data is subgroup Ic, which is located along the same line of 
sight of the ONC a nd is thought to be as young as 2 Myr old 
dBrown et al.l fl994). Thus, we expect little differences between 
the ONC and subgroup Ic isochrones, the main source of scatter 
being the smaller distance modulus, providing photometry up to 
0.3 magnitudes brighter. 

Considering that subgroup Ic is, at most, as reddened as 
the ONC by interstellar extinction, we expect that the isochrone 
shown in Fig. [TJ also applies to the CMD of subgroup Ic. Thus, 
we still expect to find in each neighborhood the two main popu- 
lations discussed above, the least reddened containing both ONC 
and subgroup Ic members, together with foreground field stars. 
Therefore, the presence of spurious sources belonging to the Ic 
subgroup should not affect our computation of the OMC-1 ex- 
tinction map. 

3.3.2. Contamination by extragalactic sources 

In deriving the OMC-1 extinction map, we model the back- 
ground sample with the reddened galactic model but neglect any 
contamination by extragalactic sources. We now address the re- 
liability of this assumption. 

To compute the extragalactic counts model reddened by the 
OMC-1, we use a deep NIR catalog of galaxies based on SOFI 
observations in an area of 340 arcmin 2 , centered on the coordi- 
nates RA=3 A 32'28" and DEC=-27°48'27" (kindly provided by 
T. Dahlen, private communication). At these coordinates, the in- 
terstellar galactic extinction is particularly low (SFD98 extinc- 
tion map provides Ay ~ 0.02), entirely negligible compared 
to the extinction provided by the OMC-1. Thus, assuming that 
the distribution of galaxies beyond the OMC-1 is well repro- 
duced by the SOFI catalog, we randomly and uniformly spread 
the sample of galaxies over the field. Then, assuming that the 
main source of extinction is provided by the OMC-1, we add the 
proper amount of extinction based on our map, converting from 
vi sible to the infrared w avelengths via the extinction curve given 
bv lCardelli etaD dl989h : 

Aj = 0.288Ay; A H = 0.182Ay; A Ks = 0.118Ay (1) 

We repeat this procedure 1,000 times, and we average the 
sample of 1,000 extragalactic CMDs, in order to minimize sta- 
tistical fluctuations. We find that the reddened extragalactic pop- 
ulation accounts for just ~ 1 % of the observed background sam- 
ple in Fig.[TJdow n to H~17. This is roughly the same fraction of 
galaxies found bv lRicci et alJ d2008l) in their analysis of the deep 
ACS images of the Orion Nebula, and confirms our hypothesis 
of negligible extragalactic contamination. We also remark that 
we removed all the galaxies resolved by ACS from the analyzed 
sample, thus further minimizing the incidence of extragalactic 
sources. 
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Fig. 4. Left panel - Extinction map of the OMC-1 derived as described in the text. The color bar shows the correspondence between 
the gray scale and visual extinctions. We also mark the boundaries of the three regions outlined by RIO. Central panel - ISPI K$ 
mosaic of the surveyed ONC. The overplotted contours represent the Ay=10 (yellow line) and the Ay=20 (red line) levels. Right 
panel - The error map for our derived OMC-1 extinction map. The OMC-1 extinction map and the corresponding error map are also 
available in electronic form at the CDS . 



4. The Orion Nebula extinction map 

It is known that most of the extinction toward the Trapezium 
cluster arises at the interface between the M42 HII region and 
the neutral diffuse matte r (seelO'Dell et al.ll2009L and r eferences 
therein). As indicated bv lO'Dell & Yusef-Zadehld2000l) . internal 
extinction may also play a role. Aiming at mapping the extinc- 
tion provided by the diffuse matter in front of the OMC-1 (the 
"Orion Nebula extinction" hereafter), we focus now our atten- 
tion on the sample of stars which are candidate members of the 
ONC. All these sources have been detected in the J, H andXj - 
bands. 

We face now a different set of problems. Whereas one can 
safely assume that the 2 Myr isochrone provides a relatively ac- 
curate locus for the dereddened stellar photospheres, one has to 
consider the spurious presence of circumstellar matter, disks in 
particular, which may affect the brightness (and to a lesser extent 
the IR colors) of the sources. Depending on the typical disk flar- 
ing angles, which are enhanced b y the UV radiation within the 
HII region dRobberto et al.l I2002T) the sample of edge-on disks 
may be relevant. We shall assume that they still represent a small 
contamination, since they tend to move the sources beyond our 
Ks = 15 limit for cluster membership. Probably more important is 
the fact, also found by RIO, that ~1 1% of the point-like sources 
in the ISPI catalog do not show NIR colors compatible with red- 
dened photospheric colors. Their position in the (J-H,H-Ks) di- 
ag ram is compatible with the reddened CTTSs locus introduced 
bv lMever et al.l(ll997l) (Fig.|5j, i.e. their photometry is indicative 
of strong exc ess emission from circumstellar disks. Following 
Mever (1996), we split our candidate members sample in two 
subgroups. If the observed colors are compatible with reddened 
CTTSs colors, then we deredden the observed colors taking the 
CTTSs locus as reference. Otherwise, in order to minimize the 
effects of eventual NIR excesses from circumstellar disks, we 
deredden our sources in the {JJ-H) diagram taking the 2 Myr old 
isochrone shown by RIO as the reference locus. In both cases, 
our dereddening algorithm moves the observed photometry to 
the reference locus (either the CTTSs locus or t he isochrone) 
along the reddening direction, computed using the lCardelli et al.l 



(1989) reddening law. Figure[5]graphically shows our deredden- 
ing algorithm. 

The set of Ay estimates is then locally averaged following 
the same strategy described in Sect. 13.21 In this case, we adopt 
a maximum angular resolution of 1,500 pixels (~7.5') to allow 
our algorithm to find a good number of cluster members in the 
outskirts of the Nebula, where the surface density of stars is gen- 
erally low (see Fig. The resulting extinction map, with the 
associated error map, are shown in Fig. [6] 

We find that in the direction of the Trapezium cluster the ON 
extinction is generally Ay <3, with a clear increase (Ay >6) in 
correspondence of the dark ridge along the north-east edge of the 
HII region, which is therefore a foreground structure. The ridge 
is part of a larger bow-shaped feature extending from the east to 
the north of the ON. These findings are fully c onsistent with the 
map derived by lO'Dell & Yusef-Zadehl (l2000h combining radio 
and optical measurements over an area 360"x425" wide around 
e'Ori-C. 



5. Discussion 

5.1. Structure of the OMC-1 beyond the Orion Nebula 
Cluster 

Alongside with the previous studies illustrated in Sect. |2l our 
map (Fig. lU shows that the OMC-1 surface density generally 
increases towards the Trapezium. In particular, we find that the 
Trapezium cluster is located in front of a high-extinction region 
(Ay > 30) extending ~10' to the north of fl'Ori-C. This region 
is delimited to the south-east side by a sharp edge: the optical 
thickness of the OMC-1 decreases steeply by ~20 magnitudes 
in a few arcminutes. This edge corresponds to the Orion Bar, 
the bright feature directly discernible in the ISPI image (right 
panel in Fig. HJ. The steep drop of extinction can be confirmed 
even by a direct inspection of the images. It is quite evident that 
the density of faint and red stars suddenly increases to the south 
of the Orion Bar. A similar extinction drop is found also to the 
north-eastern edge of the extinction peak, in correspondence of 
the dark structure, known as the Dark Bay or Fish Mouth, seen 
in absorption in the optical images of the ON. Elsewhere, the 



6 



G. Scandariato et al.: The extinction map of the OMC-1 molecular cloud behind the Orion Nebula 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

H-K s J-H 

Fig. 5. Left panel - Example of the dere ddening algor i thm, coded with the black dashed line, for a star with colors compatible with 
the reddened CTTSs locus provided by M ever et all (Il997l) (the green solid line). The blue solid line is the 2 Myr old isochrone 
suitable for the ONC. The red dashed lines are parallel to the reddening direction and bracket the regions occupied by the previous 
loci reddened by any Ay. Right panel - Example of the reddening algorithm for a star with NIR colors compatible with reddened 
photospheric colors. The blue line represents the isochrone in the (JJ-H) diagram. 




Fig. 6. Left panel - Map of the average extinction affecting the ONC candidate members. The color bar shows the correspondence 
between the gray scale and visual extinctions. We also mark the boundaries of the three regions outlined by R10. Central panel - 
ISPI Ks mosaic of the surveyed ONC. The overplotted contours represent the Ay-3 (yellow line) and the Ay-6 (red line) levels. 
Right panel - Uncertainty on our derived ON extinction map. The Orion Nebula extinction map and the corresponding error map are 
also available in electronic form at the CDS. 



extinction smoothly decreases down to Ay > 4 with increasing 
distance from the Trapezium cluster. 

On larger scales, the OMC-1 extinction map shows a north- 
south pattern. The elongated extinction ridge is distributed over 
the full extent of the survey, reachi ng the OMC-2/3 star form - 
ing region to the north of the ONC dPeterson & Megeathl l2008). 
The extinction map thus follows t he dense filament traced by the 
molecular column density data of iGoldsmith et all dl997l) . 

To compare our extinction map with the one of SFD98, we 
degrade our spatial resolution down to their value. The pixel- 
by-pixel ratio r^ v = Ay(S FD)/Av(IS PI) is plotted in Fig.|7]as 
a function of AyilSPI). For low extinction values, the SFD98 
values generally overestimate the ext inction by a factor of 1.5 - 
2, consistently with what fo und by lArce & Goodman (1999). 
Moreover, consistently with [Dobashi et all (120051) . fee SFD98 
values are systematically larger than ours, the ratio increasing in 



the 3 < rA v < 5 range with the Av(IS PI) extinction. This is in- 
dicative of the fact that SFD98 map is not accurate either in high 
extinction regions or in regions with high extinction gradients. 



5.2. Structure of the Orion Nebula 

As proposed in the 3-D model of the ON bv lO'Dell etaP (l2009h . 
the boundary between the HII region and the bulk of the neutral 
matter is made up by a geometrically thin and optically thick 
shell swept up by stellar winds. Given this model, we argue that 
the extinction in the Dark Bay direction and the broad feature 
around it are spatially correlated, the former being a clump of 
high column density gas inside the remnant of the neutral shell 
located between th e observer and the Trapezium (see Fig. 4 in 
lO'Dell et al.ll2009h . 
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Fig. 7. Diagram of the ratio r& v as a function of Ay(ISPI) com- 
puted over our extinction map. The solid line is a binned sample 
median, aimed at stressing the overall trend of the ratio. This 
diagram shows that for low extinctions our measurements and 
the SFD98 extinction are comparable, and that the latter be- 
comes up to 3 times larger for higher extinctions, as found by 
iDobashiet al1d2005l) . 

On larger scales, the ON is never thinner than Ay ~2 and 
this is compatible w ith the presence of a foreground neutral veil 
dO'Dell et al.ll2009T) . Furthermore, a north-south pattern is evi- 
dent, similarly to what is found for the OMC-1 extinction map 
(Fig. 3]i. This pattern rou ghly follows the north-sou t h elon gation 
of the ONC observed bv lHillenbrand & Hartmannl 0998), indi- 
cating that the clust er is still partially embedd ed in the OMC-1. 
As pointed out bv lHartmann & Burkertl d2007l) . this suggest that 
it might be better to view the ONC as a dense, moderately elon- 
gated cluster embedded in a larger scale filamentary stellar distri- 
bution rather than a system having a s moothly varying ellipticity 
as a fu nction of radius. As indicated bv lHillenbrand & Hartmannl 
(fl998h . the close association between the distribution of young 
stars and the elongated diffuse structures suggests that no dy- 
namical relaxation has taken place so far, i.e. the stellar popula- 
tion has not dynamically adjusted from its initial state. 

6. Conclusions 

In this paper we have used the most recent wide-field NIR pho- 
tometric catalog to discriminate between two main families of 
stellar sources in the ON region: 1) ONC members, distributed 
over an area ~30'x40' roughly centered on O'Ori-C; 2) back- 
ground stars. The two samples have been statistically analyzed 
to derive the extinction maps of the OMC-1 (Sect. [3T2b and the 
foreground ON (Sect. |4|, with angular resolution <5' and <7.5' 
respectively. 

Our results show that the OMC-1 generally accounts for the 
largest amount of extinction, typically Ay >6, steeply rising up 
to Ay >30 in the direction of the Trapezium cluster and dramat- 
ically dropping south of the Orion Bar. 

The extinction toward the ON is lower, Ay <3, with a peak 
value Ay ~6 in the direction of the D ark Bay feature. Our find- 
ings agree with |Hillenbran dl d 1997b . iMuench et al.1 d2002l) and 
iDa Rio et alJ yOlO) who derived that the ONC members are typ- 
ically extincted by Ay <3. 

To gether with lArce & Goodman! ([1999) and iDobashi et alJ 
(2005), we find that the OMC-1 extinction map proposed by 




20= 15° 
Right Ascension 



Fig. A.l. Completeness contours forKs-ll in the neighborhood 
of the Trapezium stars. The contours represent the 10%, 30%, 
50%, 70% and 90% levels with increasing distance from the in- 
ner cluster respectively. This figure shows that the iso-sensitivity 
contours roughly follow the iso-brightness levels of the image, 
indicating that the completeness of the survey is strongly af- 
fected by the background nebular emission. 



SFD98 is overestimated by a factor of~3-5, especially in the op- 
tically thickest regions. Together with the comparability bet ween 
our ON extinction map and the results of |Q' Dell & Yusef^Z adeh 
(2000), this supports both the robustness of our statistical ap- 
proach and the general validity of our results over the full spatial 
extent of the Orion Nebula Cluster. 

Our derived maps for the OMC-1 and the Orion Nebula are 
available in electronic form as FITS files at the CDS. 



Appendix A: Completeness 

To determine the completeness of the R10 photometric catalog, 
we use the artificial star experiment performed by R10, increas- 
ing the number of tests up to 10 4 per magnitude and per field. 
Our completeness estimation turns out to be a "local" procedure 
as it provides the completeness sensitivity at any given position 
and for any given magnitude, computed within a circle with a 
radius of ~1', i.e. lower than the typical size of the nebular struc- 
ture present in the ISPI images. The number of tests allows us to 
estimate the 100% completeness with an error of a few 1%. 

According to R10, we find that the completeness (and the 
sensitivity) of our survey decreases with decreasing distance 
from the Trapezium. This is due to the combined effects of in- 
creasing crowding and nebular brightness close to the inner clus- 
ter. 

The results of our simulations are shown in Figure IA.1I re- 
porting the completeness trend for a Ks = 17 star in the inner sur- 
veyed region. As we stated above, we find that the sensitivity 
of the survey gets drastically shallower in the very inner region, 
decreasing from <100% down to <10% in a few arcminutes. 

Acknowledgements. We thank Prof. C. R. O'Dell for providing his useful com- 
ments on this paper. 
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